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Solving nonlinear Volterra 
integro-differential equation by using 
Legendre polynomial approximations 
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Abstract 


In this paper, we construct a new iterative method for solving nonlinear 
Volterra Integral Equation of the second kind, by approximating the Legendre 
polynomial basis. Error analysis is worked using property of interpolation. 
Finally, some examples are given to compare the results with some of the 
existing methods. 
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1 Introduction 


The area of orthogonal polynomials is an active research area in mathe- 
matics as well as with applications in mathematical physics, engineering, 
and computer science [6, 16]. Several numerical methods were used to 
solve integro-differential equations such as successive approximation method, 
Adomian decomposition method, Chebyshev and Taylor collocation meth- 
ods, Haar Wavelet method, Wavelet Galerkin method, monotone iterative 
technique, Tau method, Walsh series method and Bezier curves method 
[2, 3, 4, 6, 13, 22]. One of the most common set of orthogonal polynomi- 
als is the set of the Legendre polynomials Lo(x), £1(x),..., Lac(a), which are 
orthogonal on [—1,1] with respect to the weight function w(#) = 1. The 
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Legendre polynomials L,(x), for —1 < # < 1 and n > 0, are given by the 
forms [6, 7, 11, 12, 21] 


Lie = > (-1)* @ (a 2 n=0,1,..., (1) 


k=0 


where [n/2] = n/2 if n is even, otherwise "5+. To use the Legendre poly- 


nomials for our purposes, it is preferable to map this to [0,1]. Then we can 
also define them by the following recursive formula [11, 12]: Lo(x) = 1; and 
L(x) = 2x4 — 1 and for n = 1,2,... 


(n+ 1) Ln4i(a) = (2n + 1)(2a — 1)L, (x) — n£Ly_1(2). (2) 


On the other hand, the methods based on Legendre polynomials may be 
more appropriate for solving linear and nonlinear differential and Fredholm- 
Volterra integral and integro-differential-difference equations [5, 6, 7, 15, 18, 
21]. Legendre polynomials are examples of eigen functions of singular Strum 
Liouville problems and have been used extensively in the solution of the 
boundary value problems and in computational fluid dynamics [5, 20] . Sev- 
eral ways for solving nonlinear integro differential equations are exist, for 
example Ghasemi et al. [8] with homotopy perturbation method and in [9] 
with wavelet Galerkin method and in [10] with sine—cosine wavelet method, 
Zhao and Corless in [23] adopted finite difference method, Lepik and Tamme 
in [19] with Haar wavelet method. In this paper, by means of the matrix 
relations between the Legendre polynomials and their derivatives, the men- 
tioned methods above are modified and developed for solving the following 
nonlinear Volterra Integro-differential equation with variable coefficients 


fi(a)u(x) + fo(a)ul(a) = g(a) + a K(a,t, u(t))dt, (3) 


where u € X := C([0,1],R), f : [0,1] — R, K : [0,1]? x R > R, also 
is assumed K is a continuous function, and wu : [0,1] — R is an unknown 
function. We have obtained a solution expressed in the form 


M 
ula) & S~ anLn(2). (4) 
n=0 


Next sections of this paper are organized as follows. In Section 2, expansion 
of Legendre basis properties and matrix relations, and its discretization of a 
integro-differential equation are given. In Section 3, the convergence of the 
method is described. In Section 4, the efficiency of the method by solving 
some examples and comparison of the numerical solutions with some other 
existing methods, is shown. A short conclusion is given in Section 5. 
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2 Expansion of Legendre basis and method of solution 


If we define L(x) = [Lo(x) Ly(x) ... Lag(x)] and A = [ao ay ... ayy)” then 


M 
u(a) = , AnLn(x) = L(x)A. (5) 
n=0 


Simillary if we define L’(x) = [LZo(x) L{ (x) ... £4,(x)] we have 


M 
ula) & Yo dnl, («) =U (x)A, (6) 
n=0 


where ’ denotes the derivative with respect to x. By using Legendre recursive 
formula (1) for n = 0,1,2,...,M, we can also obtain the matrix form of the 


equation as follows 
L'(x) = L(x)", (7) 


where 2 has two forms different for odd and even values of M, that is, 
for odd values of M we have 


000 0 0) 0 0 
1 00 0 0 0 0 
0 3 0 0 0) 0 0 
ge |i 0 5 0 0 0 0 
03 0 7 2M —3 0 0 
1 0 0 0 2M—1 0 
and for even values of MZ 
0 0 0 0 0 0 0 
1 00 0 0 0 0 
0 3 0 0 0) 0 0 
1 0 5 0 0) 0 0 


From (6) and (7) we get 
u(x) & L(x)QTA. 


We use this method to approximate left hand side of Volterra integro- 
differential equation (3) as follows 
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fila)u(2) + fola)u'(a) © fi@)L(x)A + fo(a)L(X)Q7A (8) 
= (fi(e)L(x) + fo(x)L(x)Q*)A. 


Thus 
filw)u(@) + fo(w)u'(a) © S)A, (9) 
where S(x) = [s0(x) 51(2) ... sa¢(x)], and for i = 0,1,..., M, we define 
si(2) = fila)Li(x) + fala)LO9 (OT. 
To obtain a solution of the problem (3), for each x,t € [0, 1] we define 
K(a,t,u(t)) © L(x)K*LT (t) 
where K* = [kn], and 


fence < In(x),< K(a,t, u(t)), Lm(t) >> 
~e Ll? || Lmll? 


We use this method to approximate the right hand side of Volterra integro- 
differential equation (3) as follows 


g(x) + [ L(x)K*L" (t)dt = g(a) + L(x)K* [ L™ (t)dt. 


By using Legendre formulas 


we have 
g(x) + [ K(a,t,u(t))dt = g(x) + L(x)K*(Q)~'(L(x)" — L(0)*). (10) 


Let 
h(a) = g(x) + L(x)K*(Q)~!(L(x)* — L(0)*), (11) 


then from (9) and (11) we have 
S(x)A = h(a). (12) 


We can use a matrix method based on Legendre collocation points defined 
by 
i 
t= 7 i=0,1,...,M. (13) 
Now, by substituting the collocation points into Eq. (12) we have the follow- 
ing system 
S(x;j)A =h(a;) i=0,1,...,M. (14) 
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Thus, we use this numerical method to approximate the solutions of nonlinear 
Volterra integro-differential equation, which correspond to a system of (M+1) 
algebraic equations for (MM +1) unknown Legendre coefficients dg, @1, ..., @yz- 
Briefly, Eq. (14) in the matrix form is as follows 


SA =H, (15) 
where for i = 0,1,..., M@ 
S= [S(2o) S(21) oes S(am)|", 


and 


3 Error analysis 


We assume that u(x) is a sufficiently smooth function and Pyy(x) is the 
polynomial that interpolates u at points x;, 7 =0,1,...,M that are the roots 
of M +1 degree shifted Chebyshev polynomial in [0,1]. Then we have 


(16) 


where € € [0, 1], therefore 


dM+lu(x) To (@ — 2) (17) 


|u(a) — Py (2)| < maz| dxMti | (M + 1)! 


M+1 
d. 


If we assume that c is an upper bound for mar, then 


1 
(M+ 112241" 


|u(z) — Pu(2)| <¢ 
Theorem 3.1. Let uy(x) = Uj; L(x) where Um = [uo wi ... um 


1 
Um = (2m + yf u(x) Lm (x)da, 
0 
then, there exists a real number c! such that 


(19) 


: 1 
lula) — a (@)lle < ¢ Gey peRraT 


Proof. Suppose f : [0,1] > R be an arbitrary continuous function. We define 
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IIflle = i f(a) Pae. (20) 


Let Xj, be the space of all polynomials that their degrees are equal or less 
than M and f(x) be an arbitary function. Since Xj, is a finite dimensional 
vector space, f has an unige best approximation uz, such that 


les(2) —uaclle < lls —glla Yo € Xa. (21) 
In particular, we have 
1 dl 
u(x) — ue(2) 2 = | lu(e) — uae(«) Pde < i lu(e) — Parle) Pdr, (2) 
0 0 


where Py interpolates f. Thus 


ula) — wan =f Ceo eayiganrsa Ae (23) 
so 1 
|u(z) — um (2)Il2 < “Ors Die (24) 


4 Numerical examples 


In this section, several numerical examples are given to show the efficiency 
of our proposed method for approximating the solution of Volterra integro- 
differential equation by comparing with other methods. In all examples N 
denotes the number of itereations 


Example 4.1. Consider the following nonlinear Volterra integro-differential 
equation of the second kind with the exact solution u(x) = x? 


1 1 " 
(x—1)u'(#)+2u(x) = 3(a—)a?— 32+ 3.2008(0°)+ [ at? sin(u(t))dt. (25) 
0 
Comparison of the absolute errors between Block-Pulse functions method 
[1] and the proposed method for N = 7 is shown in Table 2. Also, Figure 
2 shows the comparison between exact and approximate solutions for N = 2 
and N = 7. 


Example 4.2. Consider the following nonlinear Volterra integro-differential 
equation of the second kind with the exact solution u(x) = x — x? 


3(a2 — 1)u(x) + 2?u' (x) = f(x) + in (a — t)u(t)dt, (26) 
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Table 1: Absolute errors for Example 4.1 


t BPFss method [1] | proposed method 
for N=7 for N=7 
0.0100 4.654 x 10-4 3.585 x 10-% 
0.3537 8.098 x 107° 1.726 x 10-7 
0.6101 6.675 x 107° 6.052 * 107" 
0.9500 3.581 * 10-* 2.138 x 10-* 


* N=7 
°- N=2 
real function 


Figure 1: Comparison between exact and approximate solutions for Example 4.1 


where 
f(x) = 3(a — 1)(@ — x”) + 27(1 — 22) — ta4 + 3(a +1))x3 — ail /2)e*. 


Comparison of absolute errors between CAS wavelet method [3] and the 
proposed method for N = 7 is shown in Table 2. Also, Figure 2 shows the 
comparison between exact and approximate solutions for N = 2 and N = 7. 


Example 4.3. Consider the following nonlinear Volterra integral equation 
of the second kind with the exact solution u(x) = In(a + 1) 


u(x) =e) + f at? (u(t))*dt, (27) 
0 
where 
= : a Le > Aiea ee x)ax(In(x-4 2 1 eae) 
P(e) = tg tye get gt a (et) 5 (1409 )a(in(w+1))?— 52°(11—5 ca 
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Table 2: Absolute errors for Example 4.2 


t CAS wavelet method [3] | proposed method 
for N=7 for N=7 
0.0100 3.37 x 10-3 4.8 x 10-4 
0.3446 4.72 x 1073 5.7 x 1075 
0.7075 5.87 x 1073 3.4 x 1074 
0.9178 3.42 x 107? 2.13 x 10-6 
1.0000 6.20 x 107? 5.8 x 1075 


* N=7 
° N=2 
real function 


Figure 2: Comparison between exact and approximate solutions for Example 4.2 


Comparison of absolute errors between DT wavelet method [4] and pro- 
posed method for N = 7 is shown in Table 3. Also, Figure 3 shows the 
comparison between exact and approximate solutions for N = 2 and N = 7. 


5 Conclusion 


In this paper, we have solved nonlinear Volterra integro-differential equation 
of the second kind by using Legendre polynomial. A considerable advantage 
of this method is to find the approximation of analytical solution that is a 
polynomial of degree up to N. An other advantage of the method is that Leg- 
endre coefficients of the solution can be found very easily by using computer 
programs. The convergence of this method has been presented by Theorem 
3.1. 
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Table 3: Absolute errors for Example 4.3 


t DT wavelet method [4] | proposed method 
for N=7 for N=7 
0.054 3.29 x 10-2 3.90 x 10-5 
0.600 1.49 x 107? 1.83: x10-* 
0.851 1.82 x 107! 2.88 x 107° 
1.000 4.71 x 1071 7.52 x 10-° 


* N=7 
° N=2 
real function 


Figure 3: Comparison between exact and approximate solutions for Example 4.3 
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